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Abstract 

We propose a robust classification algorithm for curves in 2D and 3D, under special and full 
groups of affine transformations. To each plane or spatial curve we assign a plane signature 
curve. Curves, equivalent under an affine transformation, have the same signature. The 
signatures introduced in this paper are based on integral invariants, which behave much better 
on noisy images than classically known differential invariants. The comparison with other 
types of invariants is given in the introduction. Though the integral invariants for planar 
curves were known before, the affine integral invariants for spatial curves are proposed here 
for the first time. Using the inductive variation of the moving frame method we compute 
affine invariants in terms of Euclidean invariants. We present two types of signatures, the 
global signature and the local signature. Both signatures are independent of parameterization 
(curve sampling). The global signature depends on the choice of the initial point and does 
not allow us to compare fragments of curves, and is therefore sensitive to occlusions. The 
local signature, although is slightly more sensitive to noise, is independent of the choice of the 
initial point and is not sensitive to occlusions in an image. It helps establish local equivalence 
of curves. The robustness of these invariants and signatures in their application to the problem 
of classification of noisy spatial curves extracted from a 3D object is analyzed. 

1 INTRODUCTION 

Curves and surfaces are fundamental entities in computer vision and pattern recognition. For 
example, the features of 3D or 2D objects are often spatial or planer curves, and their classification 
often reduces to a classification of curves under Euclidean, affine or projective transformations. 
A direct comparison of curves, such as shape matching, generally requires registration, and the 
ensuing complexity and difficulty in its application in many important problems, have recently 
led to a renewed research interest in transformation invariants. 

Although geometric invariants have been applied to problems in computer image recognition 
and processing for decades |21[ 122] |6j EH1 El], designing robust algorithms that are tolerant to 
noise and image occlusion remains an open problem. We start by providing a brief overview of 
various types of invariants that have appeared in computer vision literature. Euclidian differential 
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invariants, such as Euclidean curvature and torsion for space curves, are the most classical. The 
affine and projective counterparts of curvature and torsion are well known. The dependence of 
curvature and torsion on high order derivatives (up to order 3 for the Euclidean group, 6 for the 
affine group and 9 for the projective group), makes numerical approximation of these invariants 
highly sensitive to noise, and therefore impractical in computer vision applications. This has 
motivated a high interest in other types of invariants such as semi-differential, or joint invariants 
[2"8"1 12"3"1 15] and various types of integral invariants [23 13 13 [23 02] . Integral invariants of a 
curve in the latter references depend on quantities obtained by integration of various functions 
along it. Since integration reduces the effect of noise, these invariants hold a clear advantage in 
practical applications. 

While explicit expressions for integral invariants are known for plane curves in 2D, they have 
thus far remained elusive for spatial curves in 3D, primarily on account of their computational 
complexity. With an increasing availability of 3D data acquisition systems and subsequent emerg- 
ing applications, interest in 3D analysis and hence robust integral invariants for curves in 3D is 
becoming essential. 

In [8 j a hybrid integro-differential affine invariant which only uses first order derivatives along 
with integrals were computed. Although a performance improvement over classical differential 
invariants is obtained, the presence of first order derivatives still affects the performance. 

In [9] , we obtain for the first time explicit formulae of integral Euclidean and affine invariants 
for spatial curves in 3D. Hann and Hickman [13J introduced and computed these for plane curves. 
The type of integral invariants, computed in this paper, may be compared with moment invariants 
\27\ [30] . We emphasize, however, the following difference: a moment invariant corresponds a 
number to a shape, whereas an integral invariant corresponds a curve to a curve. The standard 
action of the affine group on M 3 induces an action on curves. Following the approach of [13] we 
prolong this action to certain integral expressions, called potentials, and then compute invariants 
that depend on these integral variables. A direct extension of [13] to 3D, using a Fels-Olver 
moving frame construction [7] is conceptually straightforward, but the computational complexity 
makes the problem intractable. An inductive implementation of the moving frame construction, 
proposed in [17 , dramatically simplifies the algebraic derivations, as it allows one to construct 
invariants for the entire group from invariants of its subgroups: in our case affine invariants in 
terms of Euclidean ones. 

The integral invariants defined in [13] and [9] are sensitive to parameterization, or sampling 
of the curve in the discrete case. A uniform parameterization is required for two curves to be 
compared. In order to overcome this limitation, we propose in this paper local and global 2D/3D 
signatures for the special affine and full affine group. Signatures based on integral invariants are 
defined in an analogous way as signatures based on differential and joint invariants (see [3] for 
example) . The global signature of a curve depends on the choice of its initial point and does not 
allow a comparison of its fragments. It is therefore sensitive to occlusions. The local signature, 
although slightly more computationally involved, is independent of the choice of the initial point 
and is not sensitive to the occlusion effects in the image. It allows to establish a local equivalence 
of curves being compared. 

In Section 2, after reviewing the basic facts about group actions and invariants, we define the 
notion of integral jet bundle and integral invariants. Explicit formulae for affine integral invariants 
in terms of Euclidean for curves in 2D and 3D are given in Section 3, along with their geometric 
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interpretation. In Section 4 we define a global integral signature which classifies curves with a 
given initial point up to affine transformations. We also define a local signature that is independent 
of the initial point of a curve. In Section 5 a discrete approximation of the signature construction 
is tested on curves extracted from 3D objects. The curves are given as discrete sequences of 
points, with possibly the additive noise. The experiments show that signature construction gives 
a robust method for classification of curves under affine transformations. The method can be 
easily adopted to a smaller Euclidean group. 

2 Group Action and Invariants 

In this section we review the basic terminology for the group actions and invariants, as well as 
the concept of prolonging the action to jet spaces and the notion of differential invariants. We 
then introduce the notion of integral jet space and define the corresponding prolongation of the 
action which gives rise to integral invariants. 

2.1 Definitions 

Definition 2.1 An action of a group G on a set S is a map a: G x S — > S that satisfies the 
following two properties: 

1. a(e, s) = s, Vs € S, where e is the identity of the group. 

2. a(gi,a(g 2 ,s)) = a(g 1 g 2 ,s), for all s £ S and gi,g 2 £ G. 
For g 6 G and s £ S we write a(g, s) = g • s = s. 

Definition 2.2 The orbit of a point s £ S is the set O s = {g ■ s\g £ G}. 

Definition 2.3 A function f:S—>M.is called invariant if 

f(g-s) = f(s),yg£Gandy S £S. (1) 

Invariant functions are constant along each orbit and can be used to find equivalence classes of 
objects undergoing various types of transformations. 

Let QL{n) denote a group of non-degenerate n x n matrices with real entries. Its subgroup 
of matrices with determinant 1 is denoted by SL{n). The orthogonal group is 0{n) = {A £ 
QL(n)\AA T = I}, while the special orthogonal group is SO{n) = {^4 £ 0{n)\det A = 1}. The 
semi-direct product of QL{n) and W 1 is called the affine group: A(n) = QL{n) x W 1 . Its subgroup 
SA{n) = SL{n) x W 1 is called the special affine group. The Euclidean group is £ (n) = 0{n) x W 1 . 
Its subgroup SE{n) = SO{n) x W 1 is called the special Euclidean group. 

In the paper we consider the action of the affine group A(n) and its subgroups on curves 
7(i) = (xi(t), . . . , x n (t)), t £ [0,1] in R n by a composition of a linear transformation and a 
translation, for n = 2 and n = 3: 
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where matrix A € QL{n) defines a linear transformations and vector (vi, . . .v n ) £ W 1 defines a 
translation. 



2.2 Prolongation of a group action 

Our goal is to obtain invariants that classify curves up to affine transformations. The clas- 
sical method of obtaining such invariants is to prolong the action to the set of derivatives 
{x^[ \ . . . , a4 |fc = of a sufficiently high order I 

4\ t) = ^M, ^nr (() _^W. (3) 

Definition 2.4 Functions of {x%, . . . , x n ,x^ \ i = l..n, k = that are invariant under the 

prolonged action (£?p are called differential invariants of order I. 

For the Euclidean action on curves in 3D, the two lowest order invariants are called curvature and 
torsion, and are classically known in differential geometry. Analogous invariants for the affine and 
projective groups are also known. 

As noted in the introduction, differential invariants are highly sensitive to noise. We extend 
the approach of |13] from planar curves to curves in a space of arbitrary dimension. Let ^(t) 
parametrized by t S [0, 1] be a curve. We define integral variables 



x\ 



S ai '-» anl (*) = f t x 1 (t) a ^--x n (t) a -dx i (t), (4) 







where the integrals are taken along the curve "y(t) and ct\, . . . , a n are non-negative integers, such 
that «!+• • •+ai_i+aj+i+- • -+a n / 0. We call I = a%+- ■ -+a n the order of integral variables, and 

there are totally n (^jyn — i} + 1)) OI * variables of order less or equal to I. Integration-by-parts 
formula dictates certain relations among the integral variables, the coordinates x\(t), . . .x n (t) of 
an arbitrary point on a curve j(t), t S [0, 1], and the coordinates x®, . . . x^ of the initial point 
7(0). For example 

xf' 1 '°'-'°\t) = / x 2 (t)dx 1 (t) = x 2 {t)x l {t)-xlx\- / Xl {t)dx 2 {t) = X2(t)dx 1 {t)-x 2 x 1 -x [ 2 lfi '°>-' 0] {t). 
Jo Jo 

It is not difficult to show that there are 



m=l ' 



independent integral variables of variables of order less or equal to I. A canonical choice of such 
variables is given by: 



x^ 3l, "" LiJk \t) = I x h {t) a n • • • Xjh (t) a ikdxi(t), where ji < j 2 , ■ ■ ■ < j k , <% > and i > j x . (5) 
J 

For example variable a^ 1 ' ' ''"' ' is canonical, but x^ ' 1 ' ''"' ' is not canonical. 
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Definition 2.5 Let I 1 be an N[- dimensional space of independent integral variables of order I 
and less, then the integral jet space of order I (denoted S l ) is defined to be a direct product of! 1 



and two copies ofM. n , i.e S 



T x 



The coordinates x%, . . . x n of the first copy of '. 



represent an arbitrary point on a curve 7(i), t £ [0,1], and coordinates x®, of the second 

copy o/R n represent the initial point 7(0). 

The action @ can be prolonged to the curves on jet space as follows: 
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(6) 



[ai,...,a n ] 



(t) 



xj^(t)---x n ^-(t)dx- i (t). 



It is important that the integration-by-parts relations among the integral variables are respected 
by the prolonged action, and therefore the action on the integral jet space is S l is well defined. 

Definition 2.6 A function on S l which is invariant under the prolonged action (^) is called 
integral invariant of order I. 



By introducing new variables 



Xi(t) = Xi(t) - Xi, i = 1, . . . , n 



(7) 



and making the corresponding substitution into the integrals, we reduce the problem of finding 
invariants under the action ([6| to an equivalent but simpler problem of finding invariant functions 



of variables {X±, . . . , X n , X, 



[oi,— ,a n ] 



1 . . . n } under the action of GL(n) defined by 
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l Xnit) ) 



(8) 



X 



a%,...,a 



«J(t) 



ir(i)-..ir(t) dXi(t). 



Invariants with respect to Q may be obtained from invariants with respect to ([8| by making 
substitution ([7])j^] Invariants with respect to a very general class of actions of continuous finite- 
dimensional groups on manifolds can be computed using Fels-Olver generalization [7 of Cartan's 
moving frame method (see also its algebraic reformulation [H]). The method consists of choos- 
ing a cross-section to the orbits and finding the coordinates of the projection along the orbits 



This reduction by the group of translations can be put in the context of inductive method described in 
Appendix. We feel, however, that making this step "upfront" makes the presentation more transparent. 
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of a generic point on a manifold to the cross-setion (see Appendix for more details). It can 
be, in theory, applied to find the invariants under the action described by Q for arbitrary n. 
Hann and Hickman [13] used Fels-Olver method to compute integral invariants for planar curves 
(n = 2) under affine transformations and a certain subgroup of projective transformations. The 
corresponding derivation of invariants for spatial curves (n = 3) remained, however, out of reach 
due to computational complexity (it is often the case in the computational invariant theory that 
practical computations become unfeasible as the dimension of the group increases, despite the 
availability of a theoretical method to compute them |26(|5].) In [9], we derived, for the first time, 
integral invariants under the Euclidean and affine transformations for spatial curves using an in- 
ductive variation of the moving frame method [T7] , which allowes one to construct invariants for 
the entire group in terms of invariants of its subgroups: in our case, affine invariants in terms of 
Euclidean. Explicit derivation of invariants for curves of higher in the space of higher dimension 
(n > 3) remains an open problem, which seems at present, to be of more theoretical, than of 
practical interest. 



3 Integral invariants in 2D and 3D 

In this section we present explicit formulas for integral invariants for n = 2 (plane curves) and 
n = 3 (spatial curves) under the affine action ([6]). The affine invariants are written in terms of 
the Euclidean invariants. We discuss their properties and geometric interpretation. The inductive 
derivation of these invariants is outlined in the Appendix. 

3.1 Integral Affine Invariants for Curves in 2D 

The standard affine group action on curves in M 2 : 

( ) = ( au ai2 ) ( tl ) + ( Vl ) . ^ ( au ai2 )*<>> 

V f W / V °21 °22 J V J V V 2 J V °21 °22 J 

prolongs to the action on integral variables up to the third order. 

By translating the initial point 7(0) to the origin and making the corresponding substitution 
X(t) = x(t) — x(0),Y(t) = y{t) — 1/(0) in the integrals, we reduce the problem to computing 
invariants under the action Q with n = 2. Among 12 integral variables 

x^\t) = f x(tyY(tydx(t), j^o,i + j<3 

Jo 

Y^ j \t) = I X(tyY(t) j dY(t), i^0,i + j<3 (9) 
J 

we make a canonical choice of 6 independent: yP-> Q ], y[ 2 '°] ) y[ 3 '°] ; Y^ 2,1 \ yt 1 ' 2 !, as suggested 

by formula d5j). The rest can be expressed in terms of those using integration by parts formulas, 
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as follows: 



X [0,l] 


= XY-Y^°\ 


x[°> 2 i 


= XY 2 -2Y [1 > 1 \ 




= l . X 2 Y--Y^\ 
2 2 


X [l,2] 


= X -X 2 Y 2 - Y [2 ' 1] , 


x [0,3] 


= Y 3 X-3Y^' 2 \ 


X [2,H 


= -X 3 Y- -Y [3 ' 0] . 



(10) 



This reduces the problem to finding invariants under the following GL(2)-action on M 8 . Denote 
det.A := a\\a 22 — 021021- The action is defined by the following equations: 

X = a n X + a 12 Y, Y = a 21 X + a 22 Y, 

= (det A) Y^ + ^a u a 21 X 2 + a 12 a 21 XY +^a 12 a 22 Y 2 , 

= (det A) (a 22 Y^ + ^i^ 2 ' ^ + ^ 21 2 anl 3 + ^031(011022 + oi 2 o 2 i)X 2 y 

+a2iai2022^^ 2 + ^022 2 012^ 3 

= (det A) (a n Y™ + 2 ai 2 y [1 ' 1] ) + ^n 2 a 2 il 3 + anOi 2 a 2 iX 2 y + a 2 2 a 21 XY 2 

1 3 
+ 2 0iiai 2 a 2 iy , 

Y^A = (detA)(a 2 2y[ 1 ' 2 ] + ^a 2 1 y[ 3 '°]+a2ia 2 2y [2 ' 11 ) (11) 
+ |Ona 2 i 3 X 4 + ^a 21 (2ano 22 + oi 2 a 2 i) X 3 Y 

+ ^a 2 \a 22 (2 a i2 a 2 i + 011022) X 2 Y 2 + ai 2 a2i022 2 XY 3 + ^-ai2022 3 y 4 , 
yl 2 ^ = (det A) ^(011022 + a 12 02i)y [2 ' 1] + 2 ai 2 a 22 y [1 ' 2] + onoa^ 3 ' ^ 

+ | a ii«2i X 4 + ^ a n 021(011022 + 2 oi 2 a 2 i)^ 3 Y + ^a 12 a 21 (2 a u a 22 + ai 2 a 2 i)X 2 y 2 
+o 2 2 a 2 ia 22 XY 3 + ^ai 2 2 a 22 2 y 4 , 

yM = (detA)(a 2 1 y[ 3 '°]+3a 2 2 y [1 ' 2l +3a 11 oi 2 y[ 2 ' 1 ]) 

13 1 

+ -a 3 u a 21 X 4 + a 2 n a 12 a 21 X 3 Y + - ana 2 2 a 2 iX 2 y 2 + a? 2 a 2 i Y 3 X + - a\ 2 a 22 Y 4 . 

We restrict the above action to the subgroup 50(2) of rotation matrices by setting an = 
cos (f>, 012 = — sin 0, 021 = sin cp, a 22 = cos 4>. We use the moving frame method to find invariants 
as described in the Appendix. Computationally this reduces to the substitution an = ^,a\ 2 = 
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7")«2i = -^--,0-22 = where r = \JX 2 + Y 2 in ( fllj ). The resulting non-constant expressions 
comprise a set of generating invariants for the 50(2) action: 



X SE = VX 2 + Y 2 



r, 



y[1.0] = y [l,0] 



y jM] = U Y ™X-±Y™Y-±X 2 Y 2 ), (12) 

y[ 2 '°] = Ify[ 2 I °]x + 2r[ 1 - 1 ]y- 1 x 3 r- 2 xr 3 

SE r\ 3 3 

if/] = 1 fy[l,2] X 2_ y [2,l] XF+ i y [3,0] y 2_l x 3 y3 y 

SE r 2 V 3 12 J ' 

ytv] = 1 f y[ 2 .n (x 2 - y 2 ) + 2 y^l xy + -y I 3 - ] xy - -x 2 y 4 - -x a y 2 

SE r 2 y V 7 -r -r 3 ^ ^ 

yP.o] = I f y [3,0] X 2 + 3 y[i,2] y 2 3 y [2,i] xy _ 1 X 5 V _ 3 l3y3 _ 3 Jy5 \ ^ 
bK r 2 V 4 4 4 / 

The invariants with respect to the special Euclidean group are obtained by making a substitution 
of y = y—y° and X = x — x° in the above expressions ( |l2[ ): R We note that since the denominators 
in the above formulas are invariant, the numerators are also invariant. 

We use the inductive approach, described in the Appendix, to build invariants under the 
5L(2)-action defined by Eq.(|ll[) with the condition det A = 1. The inductive method yields 



5vl(2)-invariants in terms of 5£'(2)-invariants (12): 



v [i,o] . v [i,o] _ v [i,o] 
I SA ~ 1 SE ~ 1 2 



XY 

L.rom.. 1 



= X SE Y^ ] =Y^X-^Y- l -X 2 Y 2 , 

Y SA = Y^ ] X 2 E = Y^X 2 -Y™XY + ^Y™Y 2 -^X 3 Y 3 , (13) 

v [2,l] _ v [2,l] _ Y SE v [l,2] 

X SA — I SE [1,1] SB ' 

V [3,l] _ 1 | v |3,0] 3 yj^" 1 [2,1] 3 / Y]?£ ] \ [i )2] 



By replacing (X, Y) with (x — x°,y — y°) in Eq.(13) we return to the integral jet space coordinates. 
In particular, Y^f = Y^ - \XY = f* (x - x°) dy - \{x - x°){y - y°). 

2 The notation for invariants suggests a certain correspondence between the invariants and the coordinate func- 
tions of the integral jet space, which we make clear in the Appendix. 
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The following three special affine invariants are used in the next section to solve the classifi- 
cation problem with respect to both special and full affine groups: 

v [l,o] -Jm, f 



h 
h 
h 



y[i.o] _ *XY, 



1. 



X 2 Y 2 , 



(15) 



1 



-X 3 Y 3 . 



SA 

Y& ] =Y^X--Y™Y- 
SA 2 6 

y[12] = [1,2] 2 _ [2,1] xy l y[3 ,0] Y 2 

SA ^3 12 

To obtain invariants with respect to the full affine group we need to consider the effect of 
reflections and arbitrary scaling on the above invariants. We note that the transformation x — > Xx 
and y — ► —Xy induces the transformation I\ — > — A 2 Ii, I2 — ► X 4 l2 and ^3 — > — A 6 /3. The following 
rational expressions are thus invariant with respect to the full affine group: 



tA 



tA 



h 

k 



\X 2 Y 2 



(y[i,o] 

y[i.2] 



y[2,i] xy + iy[ 3 >°] y 2 



^x 3 y 3 



(y[i,o] _ \XYf 

The first of the above invariants is equivalent to the one obtained in |13j . 



(16) 



3.2 Geometric Interpretation of Invariants for Plane Curves 



The first two integral invariants (15) readily lend themselves to a geometric interpretation. In- 



variants I\ is the signed area B between the curve segment and the secant (see Figure [TJ. Indeed, 
the term yl 1 ' ] in the invariant I\ is the signed area between the curve 7(f) (whose initial point 
is translated to the origin) and the y-axis, while is the signed area of the triangle A. Their 
difference is the area B. Since the SA{2)- action preserves areas, I\ is clearly an invariant. 




Figure 1: Geometric interpretation of the invariants I\ 

The interpretation of I2 is slightly more subtle. Using that yt 2 '°] = X 2 Y — 2X <[1,1 ^ and 
rearranging the terms we rewrite I2 as 

I 2 = _I((x 2 y 2 - 3I7' 1 ' 11 ) + (X 2 Y 2 - 3Y XM)), where X = x - x°, Y = y - y°. (17) 




Figure 2: Geometric interpretation of the invariants I2 



Further, the curve j(t) is lifted from 2D to 3D by defining z(t) = x(t)y(t) (similarly to the kernel 
idea), and (17) is rewritten as 

J 2 = -h(XYZ-3X f ZdY) + (XYZ-3Y I ZdX)) where Z = XY = (x-x°)(y-y°). (18) 
3 Jo Jo 

The geometric meaning of (XYZ - 3X JjJ ZdY) is illustrated in Figure [ij The term Jjf ZdY 
is the signed area "under" the plane curve (Y(t),Z(t)) in the YZ-plane. Thus X J ZdY is the 
signed volume C under the surface F = ^(t) x [0, X(t)] in Figure |2] Since XYZ is the signed 
volume of a rectangular prism (C+D in Figure^, then XYZ-3X J^ZdY is the signed volume of 
the rectangular prism (C+D) minus three times the volume C "under" the surface 7(t) x [0, X{t)]. 
Interchanging X and Y we obtain a similar interpretation for XY Z — 3Y f ZdX. 

3.3 Integral Affine Invariants for Curves in 3D 

The standard affine group action on curves in IR 3 : 

«ii ai2 ai3 \ / x(t) \ ( Vl \ ( ° n ° 12 ° 13 \ 

a 2 i a 22 fl23 y{t) \ + \ v 2 , det a 2 i a 22 a 23 / 0, 

031 «32 «33 / \ z(t) j \ v 3 / \ a 3 i a 32 a 33 / 

prolongs to the action to integral variables up to second order. We translate the initial point 
7(0) to the origin, and make the corresponding substitution X{t) = x{t) — x(0),Y(t) = y{t) — 
y(0), Z{t) = z(t) — z(0) in the integrals. This reduces the problem to computing SX(3)-invariants 
under the action Q with n = 3. Among 21 integral variables 

X^ k \t) = f X(tyY(t) j Z(t) k dX(t), j + k^0,i + j + k<2, 
Jo 

y[ij,fc]( t ) = [ X(tyY(t) j Z(t) k dY(t), i + k^0,i + j + k<2, (19) 
Jo 

Z ihiM{t) = j X(tyY(t) j Z(t) k dZ(t), i + j^0,i + j + k<2, 
Jo 
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we choose 11 independent: Xl 1 - 1 ' ! , X^ , , , , , , Z^ , Z^ , Z^ , Z 

The rest can be expressed in terms of those using the integration-by-parts formula. Using the 
inductive approach, we first compute the invariants with respect to rotations SO(3). We find 
the following 8 independent invariants. We obtain S*£'(3)-invariants by replacing (X, Y, Z) with 
(x — x°,y — y°, z — z ) ( See Appendix for details of the derivation.) 



X SE = v 7 * 2 + Y 2 + Z\ 



7 [o,i,o] 

Z SE 



XYZ - 2 XZ^ 1 ^ + 2 YZ^ ^ - 2 ZY^ ^ 
2VX 2 + Y 2 + = P ' 



7 [0,2,0] v [l,0,0] , „ 7 [0,1,1] 7 [1,0,0] 

y [1,0,0] = Z R Y R +2Z R Z R ^ (2Q) 

/ z g,2,0]2 +4z p,l,l] 2 

7 [0,2,0] 7 [0,1,1] 7 [1,0,1] 7 [0,1,1] ^[0,2,0] y [0,2,0] , , 7 [0,1,1] 2 7 [1,1,0] 5,(0,2,0] 2 [1,0,1] 
[1,0,1] _ ~ ZZ R Z R Z R ~ Z R Z fl A R + 4Z R Z R ~ Z R y fl 



Z 



(7(0,2,0] (,[0,1,1] 5,(1,1,0] 7 [0,2,0] 2 5,(1,0,1] „ 7 [0,1,1] 2 v [0,2,0] „ 5,(0,1,1] 5,(0,2,0] V [1,0,1] 
[1,0,1] _ l/j R Z R L R ~ L R A R + /Z R A fl + 2Z R L R Y R 

Z [0,2,0] 2 +4Z [0,1,1] 2 

5,(0,2,0] 7 [0,1,1] 7 [1,0,1] 7 [0,1,1] 7 [0,2,0] v [0,2,0] , 7 [0,1,1] 2 [1,0,1] 7 [0,2,0] 2 7 [1,1,0] 
,[1,1,0] _ -^ Z fl Z iJ Z R + Z R Z R A R ~ 4Z ,R y iJ + Z R Z ii 



vL- 1 -,- 1 -,^ _ 
6 SE - 



Z^ 2fl]2 + 4Z [ ^ 1]2 

where expressions Z [ ^ 0] , Z [ ^ 0] , Y R lfi > 0] , Z^> 1] , Z [ ^ 0] , zg<°< 1] , Z^ , Y^ , X^ , X [ ^ 1] , X^ 
are provided at the end of the Appendix. 

We use them to construct the following invariants with respect to the SA(3) action. 

XsA = Zg'g^XsE, 

9V [1,0,1] 7(0,1,0] 97 [0,1,0] 7 [1,1,0] o 7 [0,2,0] 7 [l,0,0] 1 
V [1,0,1] _ lY SE Z SE ~ lZj SE A SE + 6Zj SE A SE 1 f 91 \ 

SA ~ " 2Z^ 2' 

5,(1,0,1] 7 [0,2,0] 2 

^[1,0,1] _ z 5£ z .sg 
SA „fn.i.nl 3 



SA 7 [o,i,o] 3 

£ SE 

We introduce a simpler notation for the special afline invariants which will subsequently be used 

3 The canonical choice dictated by § is Y^ , y' 2,0,0 ' , , Y^ , , , Z^ , Z^ , Z^ , Z^ , 

We made a computation with an equivalent but non-canonical set of variables. 
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to solve the classification problem with respect to both the special and the full affine groups: 



Ji = X SA = nxX + n 2 Z - n 3 Y, 

J 2 = -4 (Y^f' 1] + X SA = 2 n 2 (XYZ 2 - 3 Z^X + 3 YZ^ - ZZ^ 1 ' - 2 ZY^) 

+ n 3 (2XY 2 Z + 3XZ^-3ZX^-4YZ^-2YY^) (22) 

- 2n 1 (3ril 1 ' ' 1 l-3Zl[ 1 - 1 - ] + JZ' 1 ' 1 '°]-iy' 1 ' - 1 l) 



T _ 27 [1,0,1] 



3 



where X = x- x°,Y = y - y°, Z = z - z° and m = ^ - Z^\n 2 = - Y^°\ 



,n 3 



XA — ^[1>0,0]_ eX p ress i on f third invariant in terms of the original integral variables is 
too long to be included. 

To obtain the invariants with respect to the full affine group we need to consider the effect 
of reflection and scaling on these invariants. For A 6 1 scaling (x,y,z) — > (Ax, Ay, — \z) induces 
scaling J\ — > — A 3 Ji, J 2 — > \®J 2 and J3 — > X e Js- We therefore obtain the following two invariants 
with respect to the full affine group of transformations: 

^4a,dj/ = | (23) 
1 J i 



3.4 Geometric Interpretation of Invariants for Spatial Curves 

The first invariant J\ may be viewed as an extension of the 2D invariant 1\ to 3D. Indeed, m, 
n 2 , and n 3 represent exactly the same area as the 2D invariant I\ (in FigurdTJ in three coordinate 
planes. They are extended from 2D area to 3D volume by multiplying by X, Z, and Y respectively. 
For example, n\X is the volume C under surface F in Figure [2] and n 2 Z and are similar 
volumes obtained by relabelling of X, Y, Z axis. Therefore, the invariant J\ is the summation of 
two volumes ri\X and n 2 Z minus the volume n^Y . The geometric interpretation of the invariants 
J 2 and J3, however, remains at the present time unclear to us. 



4 Curve Classification via Integral Signatures 



The integral invariants derived in the previous section depend on the choice of the initial point 
and the parameterization of a curve. For instance, consider a planar curve 7(t) = (1/2 sint — 



cost + l,sin t + cost — 1), shown in Figure pJ-a. A curve j(t) (Figure 

/ 2 1 \ 

by a special affine transformation ( 1 K ) . A curve J a (t) (Figure 



by a full affine transformation 



2 1.5 
2 
3 



b) is obtained from 7(t) 

c) is obtained from 7(t) 



The integral invariants I\ and I 2 for curves in Fig. [3]-a and Fig. [3}b with a matching parame- 
terization coincide and are shown in Figure [4]- a and Figure [4]-b. 

As illustrated in Figure |4]-c and Figure (4j-d, these invariants change under reparametrization 
r = y/t + 1. Therefore the graph of invariants with respect to an arbitrary parameter can not be 
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Figure 3: 




(a) curve j(t) (b) special affine transformation j(t) full affine transformation 7 a (i) 



used for curves comparison. In theory one can achieve a uniform affine invariant curve parame- 
terization by using an affine analog of the Euclidean arc-length parameter da = K 1//3 (is, where k 
is Euclidean curvature and ds is Euclidean arc-length. We would like however to keep our meth- 
ods derivative free. Even when the uniform parameterization is achieved, the dependence of the 
invariants on the choice of the initial point presents another comparison challenge for matching 
closed curves, or for matching parts of the contours. 

The signature construction, proposed in this section, leads to classification methods which 
are independent of parameterization and of the initial point. Inspired by signatures based on 
differential invariants [3], we use integral invariants to construct two types of signatures that 
classify curves under affine transformation: the global signature and the local signature. Global 
integral signature is independent of parametrization, but is dependent on the choice of the initial 
point and can not be used to compare partial contours. Local integral signature is independent of 
both the initial point and parametrization. They can be used to compare parts of the contours and 
therefore can be used on images with occlusions. As our experiments illustrate they are slightly 
more sensitive to noise than global signatures, but still provide robust classification results. 



4.1 Global Integral Affine Signature 

A global integral signature of a curve is the variation of one independent integral invariant, 
evaluated on the curve, relative to another. If a curve is mapped to another curve by a group 
transformation, their signatures coincide independently of the selected parametrization. The 
global signature, however, does depend on a choice of the initial point. 



4.1.1 Global affine signature for curves in 2D 



The special affine signature of a plane curve j(t) is constructed by, first, evaluating invariants I\ 
and I2 in Eq.(13) on this curve, and then plotting the parameterized curve (I\(t), I^t)) m F° r 
instance, the signature of the planar curve j(t) shown on Figure [3ja is a plane curve in Figure [5j 
Moreover, the signature of the curve 7(i), related to j(t) by an affine transformation (Figure K3[ 



as well as their reparametrization 7(1") and 7(1") coincide with the signature of j(t). 

Similarly, a full affine signature can be defined as a parameterized plane curve (l£,l£) defined 
by invariants in Eq.(16). Alternatively we use the two special affine invariants to cancel the effects 
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1 2 3 4 5 6 7 

Parameter t 



(a) Invariant I\ for j(t) and j(t) 




1 2 3 4 5 6 7 

Parameter t 



(c) Invariant I\ for 7(t) and j(t) 



1 2 3 4 5 6 7 

Paramter t 

(b) Invariant I2 for j(t) and 7(t) 




1 2 3 4 5 6 7 



Paramter t 

(d) Invariant I2 for 7(1") and 7(f) 



Figure 4: Dependence of invariants on re-prametrization: r = \/t + 1 



of reflections and arbitrary scalings, 

ii(*) = J ^t T , = ^tW- (24) 

max t |Ii| maxi(/^) 

Both invariants are reduced relative to the range of 1 1\ \ . It is not difficult to show that maxt \ =0 
on 7 if and only if 7 is a straight line. In this case the affine signature does not exist, but straight 
line regions can be easily detected by other means. The range of I\ is from to 1. The full affine 
signature of a plane curve 7(f) is obtained by, first, evaluating 1\ and I2 on this curve and then by 

plotting the parameterized curve (^h(t), hit^j in R 2 . For example, curves in Fig. ji|a and Fig.|i|c 
are related by a non-area-preserving affine transformation. Their full affine signatures coincide as 
shown on Figure [6} 



4.1.2 Global Affine Signatures for Curves in 3D 

To construct special affine signatures for spatial curve we use invariants J\ and J2 given by 
Eq.(22). Similarly to 2D case, the special affine signature of a spatial curve (3{t) is obtained by, 
first, evaluating J\ and J2 on this curve, and then plotting the parameterized curve (Ji(t), J2(t)) 
in R 2 . 



14 



-2.5 -2 -1.5 -1 -0.5 0.5 

Invariant I 



Figure 5: Signatures for "f(t), 7(t),7(r), and 7(7") coincide 

0.7 I 1 ! ! I 1 




0.2 0.4 0.6 0.8 1 



Reduced Invariant I 

1 



Figure 6: Full affine signatures for curves 7(f) and j a (t) coincide 



For example, the signature of a spatial curve f3{£) = (sin t — 1/5 cos 2 £ + 1/5,1/2 sin t — cos t + 
l,sin 2 i + cost — 1), shown in Figure [TJ-a, is the plane curve shown on Fig ^ A curve /3(i) is 

' 0.3816 0.7631 1.1447 ' 
obtained from f5 by a special affme transformation I 1.9079 1.5263 2.2894 



2.6710 3.0526 3.4341 



is obtained from (3 by a full affine transformation 





2 3 \ 






(i 


5 6 


. As Fig 


8 




8 7/ 







A curve /3 a {t) 
illustrates, the special 



affine signatures of j3(t) and f3(t) coincide. 

Similar to the 2D case, the full affine signature for curves in 3D is obtained by reducing special 
affine invariants J\ and J2 by the range of \J\\. 



Ji(i) 



\Mt)\ 
max[ I J\ I 



J 2 (t) 



max t (Jf) 



(25) 



The full affine signature of a spatial curve (3{t) is obtained by first evaluating J\ and J2 on 
this curve, and by then plotting the parameterized curve (jl(t), J 2 (*)) in M 2 . The full affine 
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Figure 7: 




(a) original curve (3{t) (b) special affine transformation (3{t) (c) full affine transformation (3 a (t) 




-6 -5 -4 -3 -2-1 1 



Invarinat J. 



Figure 8: Signatures for (3(t) and f3(t) coincide. 



signatures of /?, f3 and j3 a coincide as shown in Fig [9} 

The advantage of global signatures is their independence of parametrization, whereas the 
result of evaluation of invariants J\ and J2 on a curve depends on the choice of parametrizations 
similarly to I\ and I2 in 2D case. The disadvantage of global signatures is in their dependence 
on the choice of the initial point of a curve. The local signature construction in the next section 
overcomes this dependence. 



4.2 Local Integral Affine Signatures 

The signatures defined in the previous section can not be used for the classification unless the 
initial point of a curve is known. This becomes an obstacle for comparing closed curves or for 
matching parts of contours. For illustration, let us choose two different initial points, black circle 
or red star, on the planar curve in Figure [To} The resulting global affine signatures are different as 
illustrated in Figures 11 a and Figures [TT]-b. We overcome the dependence on the initial point by 



introducing local signatures. To proceed with the construction of the local signature, we replace 
the integration from the initial point with integration on local segments. To retain affine invariant 
properties of the signatures, we need to partition the curve in an invariant manner. Such partition 
can be achieved using the notion of affine arc-length from classical differential geometry. Our goal, 
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Reduced Invariant J 



Figure 9: Full affine signatures for (3(i), f3(t), (5 a {t) 



0.5 


-0.5 



X 

Figure 10: A planar curve with two different choices of the initial points 

however, is to propose a derivative free method, and so we use the lowest order integral invariants, 
namely I\ for plane curves and J\ for spatial curve to obtain an equi-affine partition of a given 
curve. The details are described in the following subsections. 

4.2.1 Local Affine Signatures for Curves in 2D 

We will use I\ to partition a given curve into equi-affine sub-segments. Assume that 7 is 
parametrized by t G [0, 1]. For this purpose we define an evaluation of invariants on sub-segments 
of 7. Recall that the integration in the integral variables is performed from the initial point 7(0) 
to a current point on the curve 7O). For instance, h(t) = £XdY -\XY, where X = x(t)-x(0) 
and Y = y(t) — y(0). Thus Ii(t) is a function from [0, 1] to K. 

We define the evaluation of an invariant on sub-segments of 7 by treating the starting point 
of a segment as the initial point, and computing the value of integral variables at the end point. 
In particular, for a sub-segment defined by the parameter range [p, q] C [0, 1], we may compute 
the localization ij p '^ = f£ (x(t) — x(p)) dy(t) — ^(x(q) — x(p)) (y(q) —y(p)), and similarly for 
invariants I2 and I3 defined by Eq. (|13|) . We note that the evaluation of an invariant on a segment 
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Invariant I 



-1 -0.8 -0.6 -0.4 -0.2 

Invariant I 



(a) Initial point is the black circle in Fig, 10 (b) Initial point is the red star in Fig 



10 



Figure 11: Global signatures for the same curve with two different choices of the initial point 



is a real number. 

We choose a sufficiently small A > and define an equi-affine partition = to < t\ < 
ijv = 1 of the curve j(t), t £ [0, 1] into sub-segments by the condition 



< 



liti] 



A. 



In practice we choose A proportionally to the maximum of the absolute value of Ji, i.e., we 
choose an integer M and set A = ma ^l J *l . Note that the total number N of segments that we 
obtain in general differs from M. The local discrete special affine signature of 7 is defined by an 



evaluation of I2 and ^3 on the intervals 
rf*- 1 '**^ i = 1..JV. Figure 



r[ti-l,tj] 
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ti-i,ti], i = 1..N, that is a set of points with coordinates 
illustrates that the local discrete special affine signature 




Local Invariant I 



Figure 12: Local special affine signature for the curve shown on Figure 10 



for the curve shown in Fig [TUJ with different starting points coincides. To obtain local discrete 

i = l..N. 



( K 



affine signature of 7, we use reduced invariants, that is we plot I ^^72^ 1 ^ a ^| j3 | 
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4.2.2 Local Affine Signatures for Curves in 3D 



For a spatial curve we proceed in a similar manner as for plane curves. We use invariant J\ to 
partition a curve ^(t),t £ [0,1] into ./V sub-intervals denned by a = to < t\ < ■•■ < tN = b 
such that jj' 1 - 1 '* 1 ] = A, i = 1..JV, where A > is proportional to the maximum of the absolute 
value of J\ . We define a local special affine signature by evaluation of J2 and J3 on the intervals 



1..N, that is by a set of points on the plane with coordinates ( J, 



-ljti] j[U— 
, J3 



1..2V. 

Figure 13 -b shows the local special affine signature for a curve shown on Figure 13 -a. The 
signature does not depend on our choice of initial point. 




-0.01 
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-0.05 
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-0.07 
-0.08 
-0.09 




Local Invariant J„ 



(a) A curve with two different choices of an initial point (b) Local special affine signature 
Figure 13: A curve and its local special affine signature 



5 3D Object Classification 

The features of many computer vision and pattern analysis problems are spatial curves. We 
can hence view the classification problem as that of assorting the similarity of curves in 3D, in 
particular, when subjected to affine transformations. In this section, we apply integral special 
affine invariants J\ and J2, the global special affine signature, and the local special affine signature 
to classify curves in 3D under special affine transformations. The performance of each of the 
proposed methods is evaluated. Applying these invariants to classification of 3D objects based 
on a set of characteristic spatial curves is in line of pQ, and will be considered in subsequent 
publications. 



5.1 Experimental Design 

The Princeton Shape Benchmark [3T] provides a repository of 3D models. A subset of three 
models are shown in Figure 14 We extract a total of 100 characteristic curves, and each of them 
are re-sampled to 5000 points with the same arch-length. We applied to each curve 9 randomly 
generated 3D special affine transformations as shown on (Figure 15). To make this problem even 
more challenging and to illustrate the noise sensitivity of the proposed approach, gaussian noise 



19 




Figure 14: 3D models from The Princeton Shape Benchmark 



with distribution iV(0, a 2 ) is added to each of the variations. We therefore obtain a classification 
set of 900 curves that has to be separated into 100 equivalent classes under affine transformations. 
The training set consists of 100 original curves without any noise and transformation. The dis- 
crimination power and sensitivity to noise are analyzed using the error rate of classification. We 
implemented a Nearest Neighbor (NN) Classifier in a Euclidean Space using Euclidean distance. 
In order to illustrate the advantages of the signature, we design two experiments. The first ex- 
periment uses a common parametrization for both the training and testing, while in the second 
experiment, we choose two different parametrizations (samplings) for the testing data. 





200 500 600 400 300 600 500 400 450 600 



Figure 15: A curve and 9 variations of it under affine transformation 



5.2 Experimental Results 



Three experiments are carried out with different noise variance, namely a = 2 (Fig. 16), a = 1, 
and a = 0.5. The error rates of the three different sigma settings with the same parametrization 
are shown in Tabled. 

In Table I, both the integral invariants and the signatures perform well as indicated by the 
error rates. For comparison, the classical differential invariants have a classification error rate 
more then 80%, which makes the differential invariants practically useless. Since the order of 
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Y X Y x 

(a) (b) 
Figure 16: (a) A spatial curve without noise (b) with Gaussian Noise iV(0,4), 



Table 1: Classification error rate with same parametrization, same initial points 



Noise variance 


Jl 


J 2 


Global signature 


Local signature 


a = 0.5 


0.0022 


0.0472 


0.06 


0.07 


a = 1 


0.04 


0.12 


0.15 


0.17 


a = 2 


0.0789 


0.2233 


0.28 


0.32 



Table 2: Classification error rate with different parametrization, same initial points 



Noise variance 


Ji 


J2 


Global signature 


Local signature 


a = 0.5 


042 


0.61 


0.06 


0.07 


a = 1 


0.48 


0.70 


0.15 


0.17 


a = 2 


0.56 


0.83 


0.28 


0.32 



Table 3: Classification error rate with different parametrization, different initial points 



Noise variance 


Ji 


J2 


Global signature 


Local signature 


a = 0.5 


0.87 


0.95 


0.95 


0.07 


a = 1 


0.91 


0.97 


0.97 


0.17 


cr = 2 


0.94 


0.98 


0.98 


0.32 
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integral variables involved in J2 is higher than J 1; as well as the explicit form of J2 is more 
complicated than J±, the performance of J2 is not as good as that of J\. The global signature is 
constructed with both J\ and J2, and the local signature is based on J2 and J3. The performance 
of these signatures is therefore slightly worse than J2. 

If the parameterizations are not the same, the plots of invariants J\ and J2 with respect to 
a parameter can not be used for a classification purpose as illustrated in table II. Even with the 
lowest noise variance, the error rates are more than 0.4 for J\ and 0.6 for J2. However, neither 
the global signature nor local signature are affected. 

If we make an arbitrary selection of the initial points, both individual invariants ( J\ and J2) 
and global signature have poor performance as shown in Table III. Only local signature may be 
used to characterize a curve. 

As a conclusion, if the training data and testing data have similar parameterization and same 
initial point, either invariants or signatures may be used. Under different parametrization, the 
global signature is the best choice. With an unknown starting point, the local signature is the 
only solution. 



In this paper, we presented explicit formulas for affine integral invariants for plane and space 
curves in terms of Euclidean invariants. Based on the invariants, we constructed signatures for 
special affine groups and full affine groups for these curves. Although we have focused here on 
a more complex case of affine transformations, the integral Euclidean invariants, presented here, 
can be used to classify curves under Euclidean transformations |10j . 

Integral invariants are functions of the parameter, and hence depend on the parameterization. 
Global integral signatures provide a classification method independent of parameterization (curve 
sampling). A global integral signature depends on the choice of the initial point. Local integral 
signatures provide a classification method independent of the choice of the initial point, they can 
hence be used on images with occlusions and for comparing fragments of contours. They are 
slightly more sensitive to noise than global signatures. 

As an experiment, a classification of characteristic curves of 3D objects, subjected to random 
affine transformation and noise, was conducted by using individual invariants and global and local 
signatures. Integral invariants allow us to perform noise tolerant classification method for curves 
with respect to the affine transformations (for comparison: differential invariants give 80% error 
rate) . 

7 APPENDIX: derivation of invariants 

7.1 Cross-section and moving frame map 

Building on the works |12 | I16 | ITT]. Fels and Olver [JJ generalized Cartan's normalization procedure 
jl], and proposed a general algorithm for computing invariants. The Fels-Olver algorithm relies 
on a map p: S — > G with an equivariant property: 



6 CONCLUSIONS 




(26) 
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From Theorem 4.4 in [7] , it follows that such map exists if and only if the action of G is free and, 
in addition, there exists a global cross- section, i.e a subset K, C S that intersects each orbit O s at 
a unique point. Indeed, under the above assumption the map p may be defined by the condition 
p(s) ■ s £ K.. Then p(s) ■ s = p(g • s) • (g • s) is the unique point of the intersection of O s and /C. 
From the freeness it follows that s may be "cancelled" and hence the condition (26) is satisfied. 

If G is a Lie group acting smoothly on W 1 and both S C W 1 and K, C S are smooth sub- 
manifolds, then ^"-coordinate components of the projection t(s) = p(s) ■ s: S — > /C are smooth 
invariant functions, called normalized invariants. Normalized invariants contain a maximal set of 
functionally independent invariants, and have a replacement property, which allows us to rewrite 
any invariant in terms of them by simple substitution [151 E] • 

Although, a global smooth cross-section does not always exist, a local smooth cross-sectiorj^] 
passing through every point of S may be found for every semi-regular actionj^] The freeness 
assumption can be also relaxed to a semi- regularity assumption. With these weaker assumptions 
the above method can be used to construct local invariants I14j^| 

For algebraic groups acting on algebraic varieties, a purely algebraic counterpart of the Fels- 
Olver construction was obtained in |15l fT4"] . The algebraic method can be combined with the 
inductive approach described below. In some particular examples, including the 3D example 
presented here, the computation based on the moving frame map p turns out to be more practical. 

When the group G is of relatively large dimension, computation of invariants by either a 
geometric or algebraic approach becomes challenging. In [17J two modifications of the moving 
frame method were proposed to simplify the computation by splitting it into two steps: invariants 
of a subgroup A C G are first computed, and then invariants of the entire group are constructed 
in terms of those. For the problem at hand, we use one of these modifications, called the inductive 
approach, which is applicable when a group factors into a product of two subgroups. 



7.2 Inductive approach 

Definition 7.1 A group G factors as a product of its subgroups A and B if for any g £ G there 
are a £ A and b £ B such that g = ab. 

We write G = A ■ B. If in addition A n B = e, then for each g £ G there are unique elements 
a £ A and b £ B such that g = ab. 

Theorem G = B -A, such that APiB is discrete, and G acts freely on a manifold S then Vs E S 
there exists a local cross- section JCa, containing s, invariant under the action of the subgroup B. 
From Lemma 4.7 in |17j it follows that invariants of G can be constructed from the invariants of 
A using the following method. 

Inductive method: 

1. Restrict the G-action to A. Find a local cross-section ICa C S for the action of A which is 
invariant under the action of B. 

4 A local cross-section is defined on an open subset of U C S and Vs g U intersects each connected component 
of O s H U at a unique point. 

5 An action of G is called semi-regular if all orbits have the same dimension. 

6 A function /, defined on an open subset U of S, is a local invariant if Vs £ U there exists an open neighborhood 
G 3 of e S G s.t. condition |T|) is satisfied for all g g G s . 
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2. Construct a moving frame map pa- S ^ A denned by the condition pa(s) • s G Ka, Vs G 
S, by solving the corresponding equations. Composition of coordinate functions with the 
projection t(s) = Pa( s ) • s: 5 — > /Ca are invariant with respect to the action of A 

3. Restrict the action of G to the action of its subgroup B on the invariant subset ICa and 
choose a local cross-section Kb C /Ca- 

4. Construct a moving frame map : ^A — ► B defined by the condition ps(s) • z G Kb, Vz G 
/Ca, by solving the corresponding equations. 

5. The G-moving frame map p: S — ► G is defined by /o(s) = pb(pa(s) ■ s)pa, and G-invariants 
are the coordinate components of p(s) • s = pb(pa(s) ' s) • (pa( s ) " s ) = Ab(£a( s )) ■ '-a(s)- 



7.3 Integral affine invariants for curves in 2D 

We have a product decomposition SL{2) = B ■ A, where B = j ( ^ ^i 2 

A = SO (2) is a group of rotations. The intersection BnA = {e}, and therefore we can apply the 
inductive method as follows. 



&n > > and 



1. We restrict the «SX(2)-action ( 11 ) to the subgroup 50(2) of the rotation matrices by setting 
an = cos</>, a\2 = — sine/), a2i = sin 0, aii = cos0. A subset Ka defined by conditions, 
Y = 0, X > serves as a cross-section on an open subset of the integral jet bundle. Moreover 



Ka is invariant under the restriction of (11) to subgroup B. 



x 



VX 2 ^Y 2 



VX 2 +Y 2 



is obtained by 



2. The corresponding moving frame map Pa( s ) = 

_ _\ ~~ VX 2 +Y 2 VX 2 +Y 2 ) 

with the condition X > (see (12)). The projection la - M 5 



solving the equation Y 
Ka, obtained by substitution pa into ( 11 ), produces a point whose coordinates are invariant 



under the action of 50(2). Non-constant normalized invariants are given by (12) and 
Y$e = is the remaining constant invariant. 



3. We now restrict the action (11) to the action of a subgroup B on an invariant subset Ka- 



We obtain the following transformations. 



Xa = h\\XA, 



Y 



[i.i] 



Y 



[1,2] 



J_ V [1>1] 
fen A ■ 
1,21 



Y 



[3,0] 



iy 

Hi A 

°ll r A 



+ 36f 2 yJ 1 ' 21 +3b n b 12 Y [ *' 1] . 



v[1.0] 

A 
v[2,0] 

A 



1,0] 



■[2,0] 



[2,1] 



611^ J + 26 12^ 
v [2,l] 9 6i2 V [1.2] 
r A ~ ] ~ z b 11 1 A 



[1,1] 



(27) 



A subset Kb C /Ca defined by the equations A^ = 1, Y 
on the subset of /C^, where Yj± ^ 0. 



[2,0] 



serves as a cross-section 
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4. This leads to the moving frame map Pb(s) 



i 

X A 



Y 



2.0] 
__1 

[1,1] ■ 



2YX' ± 'X A 

o x A 

ICb, defined by lb(s) = pb(s) • s, produces a point with coordinates 



The projection ib ■ K-a 



X b = 1, 



Y, 



[i.i] 



b 



X A Y, 



Y, 



.1.2] 



1,1] 
1 i 

[1.2] 



v [l,0] 

I B 

V [2,0] 



[1,0] 



y, 



B 

[3,1] 



r, 



[2,1] 



B 



Y 



v [2,i] y| 2 ' 01 Y \i,2] 

1 A fi 1 1 J /I 



[1,1] A ' 



-V- 



y [3.0] , 3 y{ 2 -° ] y[2,l] , 3 , , 
r A +" 2 ^[i.il "T 



2 v [i,i] A 



Y 



[2,0] 



i,i] 



1" 



1.2] 



(28) 



invariant under the l?-action (27) on ICa- 



5. Replacing coordinates Y A l ' J ^ with the corresponding Ygg given by ( 12 ) produces independent 
invariants (13). 



7.4 AfRne integral invariants for curves in 3D 

We have a product decomposition SL(3) = B-A, where B = 




b u > 0,6 22 ± 



611622 



and A = SO(3) is a group of rotations. The intersection B f~) A = {e} is trivial. We again follow 
the steps of the inductive method. 

1. We restrict the SL{3) to the action of 50(3) whose elements can be represented as the 
product of three rotations: 



/ 1 








\ 


( 


cos 4> 





sin (j. 


\ 


/ 


cos 9 


— sin0 


o\ 





COS 1p 


—simp 









1 









sin # 


cos 6 





\o 


sin rp 


cos tp 




V 


— sin<^ 





COS (j. 




V 








1 / 



A subset ICa, defined by conditions, Y = 0, Z = 0, Zqh = 0, X > serves as a cross-section 
on the open subset of the integral jet bundle where X 2 + Y 2 + Z 2 > 0. The cross-section 
ICa is invariant under the action of B. 

2. The corresponding moving frame map pa is obtained by solving the equation Y = 0, Z = 
0, £[0,1,1] = with the con dition X > 0. Explicitly 



cos 



sin ( 



\/X 2 +Y 2 ' 

y 



COS ' 



sm i 



Vx 2 +y 2 

VX 2 +Y 2 +Z 2 



, COS tp 



7 [0,2,0] 



/[0,2,0p 



VX 2 +Y 2 +Z 2 



, sin^ 



+4 (z ^ 

7 [0,1,1] 



[o,i,in 



(29) 



,[0,2,0P 



+4 (Z 



[0,1,1] 



we zj^' 1 and .zj^' 2 ' ' are given on the last page of the Appendix. The corresponding set of 



50(3) invariants is given by (20). 
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3. We now restrict the action 5L(3)-action to the action of a subgroup B on an invariant 
subset of K.A- We obtain the following transformations: 

Xa = &ii Xa, 

7 [o,i,o] _ 1 7 [0,1,0] 

Z A ~~ I A 

Oil 

7 [1,0,0] _ 1 7 [1,0,0] 012 7 [0,1,0] 

Z A ~~ I A + 7 T Z A > 

022 O11O22 



if '°'° ] = 611 6 22 y' 1,0 ' ^ - 613 o 22 Z?' 1 ^ + 611 o 23 Z' 1 ' ' ] + 612 623 , 



7 [0,2,0] _ 22 [0,2,0] 
^A " ^A - 



.[1,0,1] _ 1 7 [1,0,1] 



* 0ll o 2 * 



7 [1,1,0] _ 7 [1,1,0] 23 [1,0,1] 612 „[0,2,0] 
^ - ^ +^A + ^A , 



V [1,0,1] _ V [1,0,1] , 023 7 [1,0,1] 012 „[0,2,0] 
^ " ^ i>22 ~2bu A 



A subset /Cb C /Ca defined by equations 

^[0,1,0] = 1 ^[1,0,0] = 1 y [i,o,o] = 1 ^[0,2,0] = 1 ^[1,1,0] 

is a cross-section on an open subset of Ka- 

4. The corresponding moving frame map p_g is 

7 [o,2,o] 7 [i,o,0] 7 [o,i,0] 
, _ 7 [0,1,0] , _ Z A Z A ~ Z A 
611 - Z A > 012- -jp^gj , 

Z A 

7 [0,1,0] 7 [0,1,0] 7 [1,1,0] 7 [0,2,0] 7 [1,0,0] 

z, Z A h = ~ A Z A + Z A Z A 

022 7 [0,2,0]' " 23 7 [0,2,0] ^[1,0,1] 

A Z A Z A 



013 



7 [0,1,0] 2 7 [0,2,0] 7 [1,0,0] 7 [0,1,0] 3 7 [1,1,0] 7 [0,2,0] 7 [1,0,1] 7 [0,1,0] 2 [1,0,0] 7 [1,0,1] 7 [0,2,0] 2 
-Zj a a a Zj a + Zj a Zj a - Zj a a a Zj a 1 A + l a Zj a 

^[0,2,0] ^,0,1] ^[0,1,0] 2 



The coordinate components of the projection • s: /C^ — > /Cb 

xb = 4 ' 1,01 ^ 

7 [1,0,1] 7 [0,2,0] 2 
^[1,0,1] = Z A Z A 

B 7 [0,1,0] 3 ' 

Z A 

9V [1,0,1] 7 [0,1,0] 9 7 [0,1,0] 7 [1,1,0] o 7 [0,2,0] 7 [i,o,o] 1 
y[l,0,l] _ 2y A Z A ~ 2Z A Z A + ciZ A Z A _ £ 

B 2Z jo,i,o] 2' 

are invariant under the action of B on /C^. 
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[i,j,k] 7 [i,j,k] 



5. Replacing coordinates Y\ ,Z 



with the corresponding Y^'^ , Zg'^ 



given by (20) 



produces independent invariants (21). 
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The following auxiliary expressions were used in the paper where r = \J X 2 + Y 2 , R = \JX 2 + Y 2 + Z 2 : 
z [o,ifi] = _J_^ XYZ _ 2XZ [o,i,o] +2YZ [i,o,o]_2ZY^ ^ 
z [i,ofi] = _l_^ x 2 z _ 2XZ [im +Y 2 Z-2YZ^) 

y£,o,o] = (yx 3 + XY 3 - 2 Y 2 Y^ '^ - 2 ZYZ^ ^ - 2 X 2 Y^' ^ + 2 ZXZ^ 1 ^ 

Z ^ = L_ ( 2 YX 3 Z 2 - 6 X 3 Z^ + 6 X 2 YZ^ - 6 X 2 ZY^ 

n 6 r R V 

-6 IZ 10 - 1 ' 11 y 2 + 3 XyZZ [0 ' 2 ' 01 + 4 Z 2 y 3 X + 6 XYZX [1 '°' 1] - 6 XZ 2 X [1 ' lfi] - 6 Y 2 ZZ^ 1 ' 1 ' ^ 

_ 6 y2 zy [l ) > l] +6 y3 Z [l,0 ) l] _ 3Z 2 yx [0,2,0]^ 

Z P,2,0] = 1 X 2 Z [0,2,0] _ 2 y 2^2 z + g X y Z [l,l,0] + 3 XZX [0,2,0] + g ^[1,0,1] _ g ^y^LLO]) 

Z [1,0,1] = _J_ (_ 4 Z 2 X 4 + g Z [1,0,1] X 3 + g X 2 yz [0,l,l] _ 2 X 2 z 2 y2 + g ^2^(1,0,1] _ g XYZZ \^\ 

n 6 r l R \ 

+6XZ [1,0,1] y2 _ 3y 2 zz [0,2,0] + 6 y3 z [0,l,l] _ ^2 y 4^ 

z [i,i,0] = _J_ / g Z 2 xyx [i,o,i] _ 6 ZXY 2 X^ - 3 ZX 2 YX^ + 6 Z 2 Y 2 Y^ 
6r z R z \ 

-4 ZX 5 y + ZY 5 X + 6 X 4 ^ 1 ' 1 ' ! - 6 Y i Z^ 1 ^ - 6 Zy^ 1 ' ' 1 ] X 2 + 6 ZX^l - 1 ' 1 ' y 2 
+3 Z 2 iyz' 0A °l - 6 ZX^ 1 - 1 ' ' - 3 X 3 Y 3 Z + 12 X 3 YX^' '^ + 12 IF 3 !' 1,0 ' 1 ' 

_ 6Z y3 z [l,0,l] +6X 3y Z [0,2,0] +6X y3 z [0,2,0] + g ZX 3 z [0,l,l] 

-6 X 3 Z 3 Y + 6 X 2 Z 2 Z^ 1 ' 1 ' ^ + 6 x 2 Z 2 y [1,0 ' 1] - 3 ZY 3 X [0 ' 2 ' 0] - 3 Z 3 y 3 x) 

y[l,0,l] = 1 (_g X 4y[l,0,l] _ 12 ^17^11,1,0] + g Z 2y2y[l,0,l] _ g y4^[l,l,0] _ g zx 2y X [0,2,0] 

R Qr 2 R 2 \ 

-6 Y 2 Z^ X 2 - 12 Zyzl 1 - ' 1 ! X 2 + 6 Y 2 Z 2 Z^ - 12 ya^yiw] - 6 y^yl 1 . . 1 ! 
-6 ZY 3 X [0 ' 2 ' 0] - 12 ZX 3 X [1,1 ' 0] - 3 Z 3 Y 3 X + 9 X 3 Y 3 Z + 4 ZX 5 y - 12 ZY 3 Z [1 '°' 1] 
+12 ZX 3 Z [0 ' 1,1] + 3 XY 3 Z [0 ' 2 > 0] + 5 Zy 5 X - 3 Z 2 XYZ [0 ' 2 ' 0] + 12 ZXZ 10 ' 1 ^ Y 2 - 6 Z 2 XyX [1 '°' 1] 

+6iy 3 i' 1 ' 'H6i 3 ri[ 1 ' - 1 i +6i 2 z 2 y' 1 - ' 1 ] +3i 3 r# 2 -°i) 

X [14,0] = _J_ /_ y5x 3 Z 2 y3 ^ _ g Zy2y[ l,0,l] _ g ZX 2 y [l,0,l] _ g Z y2 z [l,l,0] _ 3 ^5y 

fl 6ri? 2 V 

+6X 3 X [1,1,0] +3Xyzz [0,2,0] +3X 2 yx [0,2,0] + g Xyzx [l,0,l] +6Xy 2 x [l,l,0] 

+3 y 3 x[°' 2 '°l - 3 X 3 Y 3 - 6 Z 2 YZV->W + 6 ^IZ' ' 1 - 1 ') , 

x [i,o,i] = J_ /_ 2 ZX 4 g X 2 X [1,0,1] 2 y 2 X 2 Z + 2 X 2 Z 3 - 6 Xy^I 1 ' 1 ' ! - 6 ZXZ^ 
R 6rR V 

-6 ZYZ^ + 2 y 2 z 3 + Y 4 Z - 3 Z[°' 2 '°l y 2 ) , 

X [0,2,0] = 1 /_ X 4 y2 _ y4 ^2 _ 3 Xy 2 x [0,2,0] _ 3 X 2^2y2 + g ^^[1,0,1] + Q X Y Z Z^ 

n 3r 2 R V 

_ 3ZX 2 Z [0,2,0] + g X 2 yx [l,l,0] _ 3X 3 X [0,2,0] + g ^[1,1,0]^ _ 
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